Negative Magnetoresistance in Andreev Interferometers 
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We consider transport in a diffusive cross-shaped Andreev interferometer geometry, both the- 
oretically and experimentally. A strong and unexpected modulation of the conductance with the 
superconducting phase is found. In particular, a reversed resistance vs. phase difference is predicted, 
where the resistance is decreased by a phase gradient. A comparison of our quasiclassical calculation 
with experimental data shows quantitative agreement if we account for diamagnetic screening. 



Quantum mechanical interference effects can be ideally 
studied in heterostructures of normal metals and super- 
conductors. Large conductance oscillations in structures 
containing loops threaded by a magnetic flux have been 
predicted and observed % |, |, |, |, |, § ■ The origin of 
these oscillation is the phase coherence between Andreev 
coupled electron-hole pairs, which can be maintained over 
large distances of the order of /Ltm, making their observa- 
tion in microstructured devices accessible PI p^. 
Since the origin of these oscillation are superconducting 
correlations, they are usually suppressed by a magnetic 
field, i.e. these devices show a positive magnetoresis- 
tance. This trend can be reversed by means of tunnel 
junctions, in which case the conductance is strongly re- 
duced due to a suppression of the low energy density of 
states leads. Obviously, the suppression of superconduct- 
ing correlations then leads to a decreased resistance . 
Thus, one would expect that a magnetic field drives the 
resistance towards its normal state value, no matter in 
what direction this actually is. 

Here we report on a theoretical calculation and an ex- 
perimental measurement of the magnetotransport of a 
purely diffusive heterostructure, which shows a negative 
magnetoresistance due to an subtle interference effect. 
This is surprising, since generally one would expect, that 
both finite energies and a phase gradient would lead to 
a destructive interference, diminishing the proximity ef- 
fect. Thus, although the conductance shows the so-called 
reentrance effect as function of temperature, one would 
not expect a phase gradient to reverse the temperature 
effect. In the present Letter we study a geometry, in 
which a separation of the energy scales responsible for 
these two effects result in a negative magnetoresistance. 

We will first describe the theoretical approach, which 
predicts a negative magnetoresistance for small magnetic 
flux. To obtain the experimentally observed fully in- 
verted resistance oscillations, we will have to include dia- 
magnetic screening in the superconducting loop imposing 
the phase difference. This allows us to obtain excellent 



FIG. 1: Geometry of experimental setup and the theoreti- 
cal model. The left picture shows a micrograph of the ex- 
perimental layout. The resistance of the horizontal Ag wire 
is measured in 4-terminal configuration by current(I)- and 
voltage(V)-probes attached to points A and B. In the mid- 
dle two superconducting mirrors C and D made from Al are 
attached by Ag wires. The superconducting mirrors are con- 
nected to loop threaded by a magnetic flux, which allows to 
vary the phase difference between C and D. Theoretically we 
model the experimental setup by the structure shown in the 
right panel. The resistors Rs, Rn and Ri model the respec- 
tive Ag-wires of length Ls, Lm and Lj. 



agreement using parameters determined from the exper- 
iment independently. We note that this agreement is 
contrary to the claim in Ref . , in which a similar ex- 
perimental layout was studied. It was however claimed, 
that the experimental observations can not be explained 
within the 'standard theory' of the mesoscopic diffusive 
proximity effect. Here, we show that these observations 
are in perfect agreement with the 'standard' quasiclassi- 
cal theory. 

The system we study is shown in Fig. |l|. The left part 
shows the experimental setup and the right shows the 
theoretical model structure. The resistance of the verti- 
cal piece of normal metal is to be measured. Each part 
(left and right) of this vertical wire has a resistance Rn 
and a length Ljv. The reservoirs inducing superconduc- 
tivity (SI and S2) are attached to the middle of the wire 
by normal resistors of length Ls and resistance Rs ■ The 
actual measurement of the resistance is performed in a 
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4-terniinal geometry by attaching current and voltage 
probes (N1-N4) through resistors of length L/ and re- 
sistance Rj. Note, that these probes have to be included 
into the modeling in contrast to normal systems. In a 
proximity effect structure, the properties of these probes 
matter, since they modify the equilibrium properties (i.e. 
the density of states and therefore the energy- and space- 
dependent conductivity of the system) . 

The calculation of the resistance of the structure de- 
picted in Fig.£^has been performed numerically, along the 
lines of Ref. using the quasiclassical formalism JT6|. 

From the spectral Usadel equation the local energy- 
dependent conductivity is extracted via 
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Here is the normal state conductivity and qr^a are 
retarded and advanced Nambu Green's functions obeying 
the Usadel equation (see Ref. [|l^ for a definition). The 
solution of the kinetic equation yields for the measured 
4-terminal resistance 
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FIG. 2: Resistance versus phase difFerence. Right plot: resis- 
tance oscillations for temperatures between 200mK and 2K in 
steps of lOOmK. Left plot: resistance oscillations for temper- 
atures from to 200mK in steps of lOmK. The form of the 
oscillations at temperatures above 200mK are as expected, 
i.e. the resistance monotonously increases between to tt. 
Below 200mK the phase dependence changes qualitatively. In 
particular for small phase difference the magnetoresistance is 
negative. 




(2) 



Similarly we can find the supercurrent in equilibrium as 
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These expressions show, that the 4- Terminal configura- 
tion indeed measures the resistance of the vertical branch 
only. However, due to the proximity effect the spectral 
conductivity a{E, x) depends also on the properties of 
the current and voltage leads. As a consequence this 
will also lead to a strongly increased supercurrent for 
temperatures at which these branches are probed, as we 
demonstrate below. 

We now discuss the concrete theoretical results. The 
parameters, that we have chosen correspond to the ex- 
perimental setup shown in the left panel of Fig. |[ The 
dimensions are — lOOOnm, Ls — 225nm and L/ = 
800nm, where L/ is taken a little bit larger than in the 
experiment to account for the proximity effect into the 
wider ("lead") regions of the experiment. For the con- 
nectors to the superconductor we have taken the distance 
between the central point and the closest point of the su- 
perconductor. Taking the experimental value of the dif- 
fusion constant D « 80cm^/s, we find for the Thouless 
energy Em = TiD/Lj^ — 57mKA:B. The branch connect- 
ing the two superconductors directly has correspondingly 
Es = hD/^L% = 280mKA:B (we took the full length, 
since this is mainly used in literature to characterize SNS 
junctions). In the calculation we have assumed a temper- 
ature independent superconducting gap corresponding to 
A/fcs = 2K. 



Let us now turn to the phase dependent resistance of 
the vertical branch, measured in 4-point configuration 
through the leads N1-N4. The resistance oscillations for 
a large range of temperatures are shown in Fig. ^. Fol- 
lowing the temperature dependence of the resistance for 
zero phase difference, we observe the usual reentrance 
behaviour, with a maximal suppression of about 13% 
at a temperature around 200mK. At a phase difference 
of TT (or any odd multiple), the proximity effect is com- 
pletely suppressed for all temperatures. At temperatures 
above 200mK the resistance vs. phase oscillation are 
similar to that predicted previously (|, 0, D and found 
experimentally |^. However, below 200mK the oscil- 
lations change qualitatively. Whereas the resistance in- 
creases monotonically starting from zero phase difference 
for the larger temperatures, it decreases first for the lower 
temperatures. Only close to half-integer phases the resis- 
tance increases again. At half-integer flux the resistance 
is always equal to the normal state value, as it should be. 

Obviously, such a strong change of the phase- 
dependence should also affect other transport properties, 
for example the temperature- and the phase-dependence 
of the supercurrent. The theoretical results are shown in 
the inset of the left panel of Fig. ||. Note that, for sim- 
plicity, we have not included the temperature dependence 
of the superconducting gap. Therefore, the supercurrent 
is finite above the experimental critical temperature of 
« lAK. The critical current strongly increases below 
a temperature close to Es- At higher temperatures the 
critical current depends only weakly on the temperature 
(the exponential decrease ^ exp{— \/ ksT / E s) sets in 
only at much larger temperatures). The strong increase 
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FIG. 3: Left panel: Supercurrent phase relation (main) and 
critical current /critical phase (inset). The non-sinusoidal cur- 
rent phase relation is plotted for temperatures from lOOmK 
to IK in steps of lOOmK (from bottom to top). For phase 
differences around 37r/4 the increase is strongest. The inset 
shows the temperature dependence of the critical current, to- 
gether with the phase at which it is reached. Both quantities 
show a strong variation below IK. Right panel: resistance 
versus temperature for integer and half-integer external flux. 
The lines are theoretical results and the symbols correspond 
to the experimental data. For integer flux both, experiment 
and theory, show the usual reentrance behaviour. For half- 
integer flux, the resistance decreases from the normal state 
resistance below a temperature of about 400mK. At lower 
temperatures it drops even below the resistance for integer 
flux. We attribute the saturation of the experimental resis- 
tance at very low temperatures to heating effects. The best 
flt of the resistance for half-integer flux was obtained for a 
screening parameter 7 = 0.005e/Gjvi5jv. 



is related to the influence of the vertical branch. As su- 
perconducting correlation become strong in this branch, 
the pair breaking influence of normal lead gets effectively 
suppressed. Accordingly the supercurrent increases dras- 
tically. The phase-difference, for which the critical cur- 
rent is reached, shifts from O.GStt up to O.Stt in the same 
temperature interval. In the left panel of Fig. ^ we show 
the current-phase relation in an temperature interval be- 
tween lOOmK and 2K. It is clearly non-sinusoidal and 
also changes qualitatively at low temperatures. 

Summarizing these theoretical results, we have ob- 
served a rather surprising phase-dependence of transport 
properties in a multi-terminal diffusive heterostructure. 
However, the experimental results, shown in the left col- 
umn of Fig. ^, are quite different. The reason, as we will 
quantify below, is the diamagnetic screening in the su- 
perconducting loop. The applied magnetic flux induces 
a supercurrent in the loop, which screens the magnetic 
flux. Therefore, the phase difference between the super- 
conducting reservoirs is not directly given by the applied 
magnetic flux. This effect depends on the self-inductance 
of loop. We will show below, that accounting of this 
leads to a fully inverted phase dependence and a very 
good agreement of experimental results and theoretical 
calculations. 

To address the screening effect we consider a super- 
conducting loop interrupted by a weak link. For the 



system depicted in Fig. |l| the weak link is the normal 
metal structure. We therefore have to assume that the 
entire phase-drop occurs over the normal metal, which is 
equivalent to say that the critical current in the normal 
metal (induced by proximity) is smaller than the critical 
current in the superconducting ring connecting the two 
reservoirs SI and S2. The energy of a superconducting 
loop, containing a weak link reads, 

Eiip) = hl{^f + Ej{^). (4) 

Here, L is the selfinductance of the loop, which depends 
on the geometry, and Ej{(p) is the Josephson energy of 
the weak link. It has to be determined from the integra- 
tion of the current-phase relation I{(p). Fluxoid quanti- 
zation requires that ip = 27r<I>tot/^*o (we restrict ourselves 
here to (/j € [0,7r], for simplicity). The total flux $tot is 
the sum of the externally applied flux <I>x and the flux 
created by the supercurrent = LI{ip). Introducing 
LPx = 27r$a;/<I>o and minimizing (^ with respect to we 
obtain an equation, which determines the phase differ- 
ence as function of the external flux 

2e 

In the limit of vanishing selflnductance (L — > 0) the phase 
difference across the junction is equal to the externally 
applied flux. This is usually the desired result. On the 
other hand, if the left hand side of Eq. (^) is not neg- 
ligible, the phase at the junction differs from the ap- 
plied flux. If e.g. the current is I{(p) = Ic sin ip and 
{2e/h)LIc > 1, we find p « ip^/ {2eLIc/fi) < p^. The 
phase is always much smaller than the applied ffux. This 
holds as long as the applied flux is less than tt. If the 
external flux exceeds this value, other solutions become 
important and the phase jumps. The main consequence 
of these jumps is, that certain intervals of phases around 
odd multiples of tt can not be reached by modulation of 
the external flux. Below, we will fully account for the 
screening, by calculating the full current phase-relation 
for each temperature. Then the solution of Eq. (^) de- 
termines the actual phase difference. 

As discussed previously accounting for the screening 
effect renders a certain interval of phase differences un- 
observable. Thus, we expect that screening suppresses 
the resistance at half integer external flux, since the ac- 
tual phase difference differs from an odd multiple of tt. 
The screening effect is directly seen in the temperature 
dependence of the resistance for half-integer flux. This 
is depicted in the right panel Fig. ^ together with the 
zero-flux resistance. The value of the selfinductance pa- 
rameter 7 is chosen to match the experimentally achieved 
minimal resistance around lOOmK. Note, that without 
screening the resistance at half-integer flux would be al- 
ways equal to the normal state resistance. Taking screen- 
ing into account strongly suppresses the resistance below 
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perimental values Rn ~ Sfi and = 57mKfcs, we cal- 
culate for the selfinductance of the loop L « 2pH. On the 
other hand, we find for a circular loop of radius R = 2/xin 
and cross section = 20 x lOOnm^ a self-inductance 
Lioop = noR{ln{8R/W) - 7/4) = 9pH. Although the self- 
inductance of a loop with the experimental geometry is 
not easy to determine, this rough estimate agrees very 
well with the selfinductance determined from our fit. 

In conclusion we have shown, that a diffusive Andreev 
interferometer can show a negative magnetoresistance 
and a non-sinusoidal supercurrent-phase relation. This 
effect has been calculated theoretically and found exper- 
imentally. We obtained excellent agreement with exper- 
imental data, if the diamagnetic screening in the super- 
conducting loop imposing the phase difference is taken 
into account. This eventually leads to a fully inverted 
('TT-shifted') resistance oscillations. 

W. B. was supported by the Swiss NSF and the NCCR 
Nanoscience. 
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FIG. 4: Resistance oscillations including screening. We com- 
pare the theoretical predictions (right column) with the exper- 
imental results (left column) using the screening parameter, 
determined from the fit in Fig. |^. In the upper row the resis- 
tance oscillations are compared for temperatures of 750, 600, 
450, 350, 250, 200, 140 mK; in the lower row for 130, 105, 
95, 85, 70, 50, 20 mK. The agreement is very satisfactory, 
although not all details coincide. 



400mK with a maximal suppression of '^7%. At lower 
temperatures, the resistance reenters, similar as the zero- 
flux resistance, back to the normal states resistance at 
zero temperature. As can be seen the resistance for half- 
integer flux is also slightly below the zero-flux resistance, 
indicating a reversal of the resistance vs. flux oscilla- 
tions. This is in agreement with the experimental data. 
The deviations at high temperatures can be attributed 
to the neglection of the temperature dependence of the 
superconducting gap. 

To confirm, that the theoretical results are fully con- 
sistent with the experimental data, we compare the re- 
sistance oscillations for various temperatures in Fig. ^ 
using the screening parameter, determined from the fit 
in Fig. ^. The similarity to the experimental curves is 
quite striking, although not all details coincide. How- 
ever, it is clearly seen, that the screening effect can lead 
to the inverted resistance oscillations ("7r-shift") at low 
temperatures, i. e. the system has a negative magnetore- 
sistance. 

Let us finally comment on the numerical value of the 
screening parameter, 7 = 0.005e/GAri?jv- Using the ex- 
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